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Abstract 

In the context of warped extra-dimensional models with all fields propagating in the bulk, we 
address the phenomenology of a bulk scalar Higgs boson, and calculate its production cross section 
at the LHC as well as its tree- level effects on mediating flavor changing neutral currents. We 
perform the calculations based on two different approaches. First, we compute our predictions 
analytically by considering all the degrees of freedom emerging from the dimensional reduction (the 
infinite tower of Kaluza Klein modes (KK)). In the second approach, we perform our calculations 
numerically by considering only the effects caused by the first few KK modes, present in the 4- 
dimensional effective theory. In the case of a Higgs leaking far from the brane, both approaches 
give the same predictions as the effects of the heavier KK modes decouple. However, as the Higgs 
boson is pushed towards the TeV brane, the two approaches seem to be equivalent only when one 
includes heavier and heavier degrees of freedom (which do not seem to decouple). To reconcile 
these results it is necessary to introduce a type of higher derivative operator which essentially 
encodes the effects of integrating out the heavy KK modes and dresses the brane Higgs so that it 
looks just like a bulk Higgs. 

PACS numbers: ll.10.Kk, 12.60.Fr, 14.80.Ec 



a mariana.frank@concordia.ca 

b n_pour@live.concorida.ca 

c mtoharia@physics.concordia.ca 



I. INTRODUCTION 



Warped extra dimensional models have become very popular because they are able to 
address simultaneously two intriguing issues within the Standard Model (SM): the hierarchy 
problem and the mass/flavor problem. They were originally introduced to treat the first issue 
[1] in a setup where the SM fields were all localized at one boundary of the extra dimension. 
Later it was realized that by allowing fields to propagate into the bulk, different geographical 
localization of fields along the extra dimension could help explain the observed masses and 
flavor mixing among quarks and leptons [2, 3]. Flavor bounds and precision electroweak tests 
put pressure on the mass scale of new physics in these models [4], but extending the gauge 
groups and/or matter content (e.g. [5]) or by slightly modifying the spacetime warping of 
the metric (e.g [6]), it is possible to keep the new physics scale at the TeV level at the reach 
of the Large Hadron Collider (LHC). 

Electroweak symmetry breaking can still happen via a standard Higgs mechanism in 
these scenarios (although it can also be implemented as as Pseudo-Nambu-Goldstone boson 
(PNGB) [7] or described within the effective theory formalism [8, 9]). As the LHC announced 
the discovery of a light Higgs-like particle of a mass around 125 GeV [10], it becomes crucial 
to have a detailed prediction of the properties of the physical Higgs particle in these models. 
The Higgs boson itself must be located near the TeV boundary of the extra dimension in 
order to solve the hierarchy problem, and so typically it is assumed to be exactly localized on 
that boundary (brane Higgs scenario). Nevertheless, it is possible that it leaks out into the 
bulk (bulk Higgs scenario), and in doing so indirectly alleviate some of the bounds plaguing 
these models [11]. 

The calculation of the production cross section of the brane Higgs in these scenarios 
has been addressed before [12-18] but we will pay close attention to the more recent works 
of [16-18]. The towers of fermion Kaluza-Klein (KK) modes will affect significantly the 
SM prediction and in [18] it was found that the Higgs boson production rate can receive 
important corrections, either enhancing or suppressing the Standard Model prediction. The 
suppression or enhancement depends on the model parameters considered, in particular 
on the phases appearing in the different Yukawa-type operators present in the 5D action. 
Previously, the analysis of [16], in which only the first few modes were considered, gave no 
contribution to the rate from the towers of KK fermions. Finally, the analysis of [17] seems 
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to indicate that with just a few KK modes a substantial effect is obtained, but of opposite 
sign as the one predicted from summing the infinite tower [18]. 

In this work we consider the effects of allowing the Higgs boson to propagate in the bulk, 
with its profile more or less localized towards the IR brane depending on the value of the 
mass parameter (3, related to the bulk mass of the 5D Higgs field. 

To keep matters as simple as possible we will set up a model containing a single family 
of up-type 5D fermions along with a bulk Higgs scalar. Generalization to a more realistic 
scenario is straight forward but we prefer to stay as transparent as possible due to the many 
subtleties involved in the calculation. 

We first compute the contribution of the complete tower of KK fermions to the Higgs 
production cross section as well as to the tree-level shift happening between the light fermion 
mass and its Yukawa coupling (leading to flavor violating couplings when considering three 
fermion families). These calculations, as outlined in [18, 20], are analytically straightforward 
and allow us to obtain simple and compact results. We then repeat the same analysis 
numerically from the point of view of an effective theory in which only the first few KK 
fermions contribute. We show that for a bulk Higgs with a thickness of the order of inverse 
TeV scale, the results obtained are the same as the results obtained by summing the complete 
KK tower (i.e. heavier modes decouple). Moreover, these results are consistent with the 
predictions obtained in [18, 20] for the specific case of a brane localized Higgs. The two 
aproaches outlined seem to give different predictions as the bulk Higgs is continuously pushed 
towards the brane. It turns out that in order to maintain the consistency of both approaches 
we need to include in the analysis the effects of a special type of higher order operators. 
After these effects are included, we will come back and address in the discussion section the 
differences among the existing calculations in the literature and stress the importance of 
including the mentioned higher order operators in the analysis. 

Our paper is organized as follows. In Sec. II we summarize the simple 5D warped space 
model used in the calculation. In Sec. Ill we present analytical results for the Higgs flavor- 
changing effects (HI A) and production (III B) , using the full tower of KK fermions. We use 
numerical methods to calculate the effects of including just a few KK modes in Sec. IV, 
both for flavor-changing neutral currents effects (IV A) and Higgs boson production (IV B). 
We include the effect of the higher order operator in Sec. V and discuss the misalignment 
between the Higgs boson profile and its vacuum expectation value (VEV) in Sec. VI. We 
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discuss the significance of our results, compare them to previous analyses and conclude in 
Sec. VII. We leave some of the details for the Appendices A and B. 



II. THE MODEL 

We consider the simplest 5D warped extension of the SM, in which we keep the SM local 
gauge groups and just extend the space-time by one warped extra dimension. 
The spacetime metric is the usual Randall-Sundrum form [1]: 

ds 2 = ^(rj^dx^dx" - dz 2 ^j , (1) 

with the UV (IR) branes localized at z — R (z — R'). We denote the SU(2)l doublets by 
Q l (x, z) and the SU (2) l singlets by W z) where i, j are flavor indices and x represents the 
4D spacetime coordinates while z represents the extra dimension coordinate. The fermions 
are expected to propagate in the bulk [2, 3]. 

The up-sector fermion action that we consider is therefore 

Sfermion = J d^xdz^ % - (Q^V A Q, - V aQ^ &) + ^<2.Q* + 

% - [U 3 Y A V A U, - V A U 3 T%) + '^UjU, + (1^ QiHUj + h.c) , (2) 

with T>a being the covariant derivative, and we have added a Yukawa interaction with a 
Higgs field H which in principle can be either brane or bulk localized. From the 5D fermion 
mass terms one defines dimensionless parameters c Ui , c qi which are a priori quantities of 
0(1). The coefficients Y*j have inverse energy units (1/y/X) since Yukawa couplings in 5D 
are higher dimensional operators. 

After separating 5D fields into left and right chiralities we impose a mixed ansatz for 
separation of variables 

q L (x, z) = q° L (z)q° L (x) + <&(z)*\(x) + ... , (3) 
q R (x, z) = q* R {z)u' R {x) + Q R (z)^ R (x) + ... , (4) 
u L (x, z) = u° L (z) q ° L (x) + UKz^Kx) + ... , (5) 
u R {x, z) = u R (z)u° R (x) + U R (z)^ R (x) + ... , (6) 

where q Q L {x) and u° L (x) are the SM fermions and ^2 r( x ) are ^ ne heavier KK modes. In order 
to obtain a chiral spectrum, we choose boundary conditions for the fermion wavefunctions 

Qr( ), «l( ), ««(++), (7) 



so that before electroweak symmetry breaking only q° L and u° R will be massless (zero modes) 
with wavefunctions: 

ql{z) = f{c q )^^z 2 - c % (8) 
u° R (z)=f(-Cu)^^z 2 ^, (9) 



where we have defined f(c) = yj 2e and the hierarchically small parameter e = R/R' ^ 
1CT 15 . Thus, if we choose c q (—c u ) > 1/2, the zero mode wavefunctions are localized towards 
the UV brane; if c q (—c u ) < 1/2, they are localized towards the IR brane. 

In order to implement minimally the Higgs sector out of a 5D scalar we use the following 
action [22] 

3 

/ j. .j4_ / n \ 

'Higgs 



Tr\V M H\ 2 - ^Tr\H\ 2 

z l 



V uv (H)5(z -R)- V IR (H)S(z - R'), 

(10) 

where /x is the 5D mass for the Higgs boson. The boundary potentials Vuv(H) and Vir(H) 
yield boundary conditions that can accommodate electroweak symmetry breaking, so that 
one obtains a Higgs VEV with a non-trivial profile along the extra-dimension. Around that 
VEV, one should then add perturbations and obtain the spectrum of physical modes, i.e. a 
SM-like Higgs boson and a tower of KK Higgs fields. The expansion should look like 

H(x, z) = Vf} (z) + hp{z)h{x) + (11) 

and we can choose the boundary conditions such that the profile of the Higgs VEV vp(z) 
takes the simple form 

v p (z) = V{fi) z 2+ P, (12) 

where (3 = a/4 + /i 2 and 



(7 , 2(1 + /3) 

' ' ~ l/ i? 3 (l - (R'/R) 2 +W) {R'y+P' 1 ; 

where f 4 is the SM Higgs boson VEV. One should note that the wave function hp(z) of the 
light physical Higgs (lightest KK Higgs field) will have the form 

so that for a light enough Higgs boson mass both profiles hp(z) and vp(z) are aligned (i.e. 
proportional to each other). 



f R,+ 

lim / h 2 (z)f(z)dz = f(R'). (16) 

?->oo J R 



The previous bulk Higgs sector is capable of reproducing the brane Higgs limit, since the 
wavefunction of the light Higgs (and its VEV) both depend exponentially on the parameter 
(3. As this parameter is increased, the wavefunctions are pushed more and more towards 
the IR brane mimicking a perfectly localized Higgs sector. 1 Indeed, the wave function of the 
Higgs can act as a brane localizer since 

lim h 2 {z) = lim v 2 {z) = 5{z - R'), (15) 

/3^oo /3->oo 

where the Dirac delta function is defined as the limit of a sequence of functions with in- 
creasing value of (3. One can easily prove that for any wavefunction f(z) (or a product of 
wavefunctions) we have 

'•/." 

There is however an issue about localizing the whole Higgs sector towards the brane since 
we just showed that only quadratic Higgs operators will "become" brane localizers. When a 
5D action operator contains more than two (or less than two) Higgs fields, the (successful) 
localization of such operators is not guaranteed. In fact in order to ensure that the 5D 
bulk Higgs scenario correctly tends smoothly to a fully localized Higgs sector, one should 
implement a prescription enforcing a precise (3 dependence on the coefficients of all operators 
containing Higgs fields. More precisely, the coefficient Y N ((3) of an operator containing N 
Higgs fields (before electroweak symmetry breaking) should behave as 

Y N (f3)=Y 1 N xf3^, (17) 

where Y X N = Y N (1). This is the only way to ensure that we can have 
fR'+ rR' 

fi^too J R — X f ,,, 

or in other words 



r R'+ f R'+ 

lim / Y N (I3) h N (z)f(z) dz = lim / Yf ^ h N (z)f(z) dz = Yf f(R'), (18) 

P^°° JR l 3 ^ 00 JR 



lim Y N {(3) h N {z) = 5{z - R'). (19) 

In particular for 5D Yukawa type couplings this prescription implies that the 5D Yukawa 
coupling will have to carry a ^/J3 dependence in order to ensure that the brane limit Yukawa 



1 Moreover the masses of the heavier KK Higgs fields depend linearly on the fi parameter and so these fields 
will decouple from the theory for very large (3. 
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coupling is non- vanishing [21] (see also [20]). But it also means that any other 5D action 
operator containing a single Higgs field would need to carry the same \[fi dependence. On 
the other hand, 5D action operators containing 3 Higgs fields (like the operator H 2 HQU) 
would have a diverging limit for (3 large unless its action coefficient Y 3 ((3) is itself suppressed 
by 1/v^. 

The previous prescription makes it technically possible to define a localized Higgs sector 
from a 5D bulk Higgs field, but it certainly seems quite contrived to appropriately fix all 
operator coefficients such that they all can give non-zero and finite contributions when the 
Higgs is localized. A brane localized Higgs sector could seem "un-generic" or "un-natural" 
if it is to be seen as a limiting case of a bulk Higgs. More details about the complete 
prescription for operators containing Higgs fields are presented in Appendix B. 

III. HIGGS PHENOMENOLOGY: ALL KK FERMIONS 

For completeness and consistency, we present first a result previously obtained in [20], 
namely the computation of the shift between the light SM fermion mass term and its Yukawa 
coupling with the Higgs field (leading to Higgs mediated FCNC when more than one fermion 
family is considered). We then calculate the coupling between the physical Higgs and two 
gluons for the 5D bulk Higgs case. 2 

We can follow two routes to obtain our predictions. The computation of the flavor 
violating couplings of the Higgs scalar with fermions will be obtained in an approach based 
on considering first electroweak symmetry breaking and then solving the 5D equations of 
motion for the fermions (i.e. the effect of the Higgs VEV is directly taken into account in the 
equations of motion and during the dimensional reduction procedure). The alternative (and 
equivalent) approach would be to consider first the dimensional reduction (i.e. obtain the 4D 
effective theory in the gauge basis), and then consider the electroweak symmetry breaking 
in the presence of the infinite tower of KK fermions. After performing the diagonalization of 
the infinite fermion mass matrix (as well as canonical normalization of the fermion kinetic 
terms) we should recover the same results. We use the first approach in the first subsection, 
and the second approach in the computation of the Higgs coupling to gluons and also in the 
2 We follow very closely the general procedure outlined in [18] and explicitly compute the prediction for the 
bulk Higgs case and then compare it with the brane localized Higgs limit that was presented there. For the 
sake of simplicity here we assume the matter fields belong to the usual SM gauge group. 
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following sections where we will truncate the infinite mass matrix in order to consider only 
the effect of the first few KK modes. 

A. Higgs Flavor violating couplings 

After imposing electroweak symmetry breaking in the Higgs sector, the four profiles 
Ql,r{z) and ul,r(z) introduced in eqs. (3) to (6) must obey the coupled equations coming 



from the equations of motion: 

-m u q L -q' R + C9 ~^ 2 qR + (^j Vp(z)Y u u R = 0, (20) 

-m* u q R + q' L + ^J^Ql + (^j vp(z)Y u u L = 0, (21) 

-m u u L - u' R + ° u ^ 2 ur + vp(z)Y* q R = 0, (22) 

-ml u R + u' L + ————Ul + (-) v p (z)Y: q L = 0, (23) 



where the ' denotes derivative with respect to the extra coordinate z and Y u is 5D Yukawa 
coupling. 

It is simple to deduce from these equations an exact expression for the mass eigenvalue 
m u in terms of the fermion profiles [20] 

m u = R 4 £dz (^(K| 2 + \q R \ 2 ) + ^^(YuURql-Y^Rul^j , (24) 

and compare it to the expression of the fermion Yukawa coupling, i.e 

yt = R 5 ( R dz^fi(Y u URq* L + Y d *q R u* L ), (25) 

JR z 

where hp(z) is the profile of the physical Higgs field. 

With these two expressions we compute the shift (or misalignment) between the fermion 
mass m u and the Yukawa coupling y\ as 

A u = m u - v 4 yt (26) 

which becomes simply 

A" = R 4 £'dz (^(M 2 + \q R \ 2 ) - 2Y: R ^q R u^ . (27) 
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In order to proceed further, a perturbative approach is used, such that we assume that 
(Y u V/lB!) < 1 where v± is the SM Higgs VEV. Knowing the analytical form of the VEV 
profile vp(z) and using the {Y u v^R') small parameter it is possible to solve perturbatively 
the system of coupled equations (20) to (23) to any order in {Y u v&B!) (see [20] for details). 
The result for the shift in the top quark Yukawa coupling is 

A* 2m 2 t 2 + c u -c q + f3 
-ix 



mv 4 t> 4 (1 - 2c q )(l + 2c u ) 



1 1 

+ 



6 + c u - c q + 3/3 5 + 2c u + 2(3 

(28) 



5 - 2c q + 2(3 4 + c u -c q + (3_ 

where we have only included the contribution from the third term in eq. (47), as the other 
terms are subdominant for light quarks, although not necesarily for the top quark. For 
clarity we omit their analytical expression here, but the complete analytical result can be 
found in [20] and in the Appendix A of this work. The shift in the Yukawa coupling has 
some dependence on the Higgs localization parameter (3 and it is shown in Figure 1 as the 
"infinite sum" result, as the procedure we followed is equivalent to diagonalizing the infinite 
fermion mass matrix in the gauge eigenbasis. 



B. Higgs production 

In this section we follow the approach of working with the infinite fermion KK modes with 
wavefunctions in the gauge basis. This is not the physical basis after electroweak symmetry 
breaking since Yukawa couplings will introduce off-diagonal terms in the infinite fermion 
mass matrix, which should be properly diagonalized in order to obtain the physical basis. 

Since the Higgs field is not charged under QCD the main contribution to its coupling 
to gluons comes from a top quark loop, as shown in Figure 2; if the model contains many 
heavy quarks the resulting cross section for the process is gg — > h is [19] 

2 

6(s -m 2 h ), (29) 

with tq = ml/ArriQ, s being the gg invariant mass squared and Q representing the physical 
fermions with physical Yukawa couplings Yq and masses thq. The form factor is given by 

A 1/2 (r) = 3 -[r + (r-l)f(r)]r-^ (30) 
9 
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FIG. 1. (color online). The shift in the top quark Yukawa coupling as a function of the bulk Higgs 
localization parameter (5. Each line represents an effective theory containing the given amount 
of KK fermions. The lower line (blue) represents the contribution from the infinite tower of KK 
modes. Apart from the direct phenomenological impact of this result, this term also affects the 
hgg coupling, as discussed in the text. The dimensionless 5D Yukawa couplings are fixed at Y = 2 
and the KK scale is set at — = 1000 GeV (the overall effect scales as Y 2 v 2 R' 2 ). 



Here we want to figure out the contribution to the hgg coupling coming from 5D quark 
doublets and a singlets, i.e. containing the SM quarks (which includes doublets and sin- 
glets (qL,UR)), along with the associated towers of vector-like KK fermions, (Ql, Ur). The 
relevant quantity to calculate is 



where i/q is the physical Yukawa coupling of the physical Dirac fermion Q and itlq is its 
mass. As stated before, it will prove useful to work in the gauge basis, and so we represent 
the Yukawa couplings between the KK fermions Ql(x) and Ur(x) in the gauge basis as 



with 




(31) 




(32) 
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FIG. 2. Loop diagram showing the contribution of the quark Qi to the Higgs-gluon-gluon coupling. 
In the SM, the dominant contribution is through the top quark due to its large Yukawa coupling 
with the Higgs boson. In RS the heavier KK fermions contribute to the coupling with potentially 
large effects, either suppressing or enhancing the SM coupling, depending on the phases present in 
the different Yukawa-type operators present in the 5D action, and on the localization of the Higgs 
(see text for details). 

Yq l u r . Its values will be obtained by performing the overlap integral of the Higgs profile 
and the corresponding bulk fermionic wave functions, i.e. 



where we have assumed that the nontrivial Higgs VEV and the physical Higgs profile are 
perfectly aligned 3 . The Yukawa couplings between different chirality KK fermions and also 
between zero modes and heavy KK fermions are obtained and written in a similar way so 
that we can write the infinite dimensional fermion mass matrix as 



where Mq = diag(MQ 1 , Mq 2 , ...) and Mjj = diag(M[/ 1 , Mu 2 , •••) are the KK mass matrices 
for the corresponding fermion fields in the gauge basis, and we have suppressed fermion 
family indices to simplify notation. From eqs. (30) and (31), we notice that in eq. (29) the 
form factors, Ai/ 2 ~ for light fermions, and A\/2 » 1 for the much heavier KK modes and 
the top quark. Therefore, separating the contribution of the light fermions from the heavy 

3 We address the case where hp(z) ^ vp(z)/vi in Section VI. 




(33) 




(34) 
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ones we write 

light v heavy v 

where in the first (second) term the sum is only over light (heavy) fermion generations. 
Noting that 

E^Ef^™"'), (36) 

heavy ^ light ^ 

where M is the fermion mass matrix given in (34), while Y is the Yukawa matrix, we have 

c hgg = TrfYM- 1 ) + — (A/ 2 (r ) - 1). (37) 
hght m ® 

We also note that Y = and since the trace is invariant under unitary transformations, 
we can compute it in the gauge basis (so we can use the fermion mass matrix in that basis). 
Up to first order in t> 4 one finds 

TVVA/r 1 ^ - din Det(M) 1 ^ 2 f KlUh/uI.Qh/q^uA 
Tr(YM ) - — — - - v 4 ^ Ji - M - [Y QliUr Y Ul]Qri — j 

Noting that the SM masses and Yukawa couplings are also modified (shifted) as [20] 



VQ I 

m Q ^ ht v A 



we can write the total hgg coupling as 

Y u i 

\ ^ Qhi^Rj U Lj Q Ri VQ I At \ /AC\\ 

c h99 = -2.4 £ m + — | ligh /i/ 2 (TQ Iight ). (40) 

where we have used equations (37), (38) and (39). As we mentioned before, the form factor 
is negligible for the light fermion generations. Therefore neglecting the last term above, and 
using (33) we have 

where the 5D bulk physical Higgs profiles can be normalized as [22] 
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with e = R/Rl ~ 10 15 being the warp factor. The sums in eq. (40) are given by [18] 
J2 9l( z )Qr( z ') _ z ,2+c *z 2_c * L„, (z'/R) 1 ' 20 " - 1" 

i=l 



.R 4 



0(z' - ^) 



e 2c q -l 



and 



E 



ujp(z)uz>(z! 



00, 



i? 4 



- ^) 



(z'/R) 



1 



-2c„-l 



(43) 



(44) 



Substitution of these sums and of the Higgs profile in Eq. (41) and assuming 4 (3 > 2 will 
finally give the total Higgs coupling for the light fermions which is given in Appendix I. If we 
assume that c q > 1/2 and c u < —1/2, which is the case for light fermions (up-like fermion), 
the expression for Ch gg can be simplified as 



Up 



v 4 Y u Y u *R' 2 



2(1 + 0) 



(45) 



^ — * " ( 2 + / 3 + Cg - Cu )4 + 2/3- 
In the case of the top quark, we have to add the contribution due to the last term in Eq. (40), 
since A 1 / 2 (Tt op ) ~ 1. Following the notation in [20], we write the additional contribution as 

^top 

llght A/2(r Qlight ) + — — 



m Q 



m t v A 



(46) 



where the first term is given by eq. (39) multiplied by the form factor, Ay 2 and last term, 
is the result of kinetic term corrections due to the shift in Yukawa couplings, which are also 
not negligible for the heavy fermions. The shift is given by 

Ar = R 4 £dz{^(\u L \ 2 + \ qR \ 2 )). (47) 

For a complete discussion on this, we refer the reader to [20]. 

So finally for IR localized fermions with c q < 1/2 and c u > —1/2 (top-like) we have 

^top 

Al/2(T Ql M ) + 



Top 

"hgg 



VQ 
m Q 



light 
v 4 Y u Y u *R' 2 



m t V4 



+ 



4 + 2/3 ' 2(3 + 5 - 2c q + 2(3 + 5 + 2c u (3 + 4 - c q + c u \ ^ 
Following our ansatz for localizing the Higgs sector, and in order to compare with previous 
brane Higgs results, we need to replace the 5D Yukawa couplings with the dimensionless 
and /9-independent couplings 

v/2(l + /3) 



Y 



-Y 



5D 



*For a completely flat bulk Higgs, 



(2 - c q + c u + (3) 

2. For any physically acceptable model (3 > 2. 



(49) 
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FIG. 3. Contribution to Chgg/ C hg^ coming from the KK partners of the "up" quark (left panel) and 
from the full top quark sector (right panel) as a function of the bulk Higgs localization parameter 
f3. Each line represents the numerical result obtained in an effective theory containing the amount 
of KK fermions indicated. The upper line (blue) represents the contribution of the infinite tower of 
KK modes (computed in the text analytically). The dimensionless 5D Yukawas are fixed at Y = 2 
and the KK scale is set at — = 1000 GeV (the overall effect scales as Y 2 v 2 R' 2 ). 

The results obtained in this section, of the contribution of a 5D top-like quark and a 5D 
up-like quark to the hgg coupling are shown in both panels of Figure 3 as the "infinite sum" 
result. 



IV. HIGGS PHENOMENOLOGY: INDIVIDUAL KK MODES 

In this section we take a different approach and compute the effects on Higgs phenomenol- 
ogy (FCNC and production cross section) due to only the first few KK fermions in the model. 
That is, we consider a 4-dimensional effective theory which contains the SM matter content, 
augmented by a few levels of KK fields. This procedure is better fitted within the frame- 
work we work in (low cut-off effective theories), the drawback being that it is not possible 
to obtain general analytical predictions in a close form. Our strategy will be to assign some 
generic values to the parameters of the model and perform the computations numerically. In 
particular we will fix the bulk mass parameters of the 5D fermions Q and U to be c u = —0.6 
and c q = 0.6 (for an up- type quark) and c u = and c q = 0.4 (for a top quark). The value 
of the dimensionless 5D Yukawa coupling will be taken to be Y^ D = 2. 

14 



A. Higgs Flavor violating couplings 



In order to evaluate the shift in the Yukawa coupling of the SM fermion (the zero mode) 
due to the presence of a finite number of KK fermions, we can simply use Eq. (39), with the 
understanding that now the sum is finite, and so we shall sum up to the maximum number 
of KK modes chosen. We are interested in computing the top quark Yukawa shift as it 
is the most interesting for direct phenomenology, and also because it will also enter in the 
calculation of the hgg coupling. We perform the sum numerically and stop the summation at 
different maximum numbers of KK fermions. The results are shown in Figure 1 in which we 
focus on the variation of the Yukawa coupling shift with respect to the bulk Higgs localization 
parameter (3 and we compare these to the results obtained in the previous section for the 
infinite KK degrees of freedom. The main observation is that for small /3, the finite sums 
are in good agreement with the infinite sum result. On the other hand for large values of 
(3 the Yukawa shift obtained from the finite sums becomes more and more irrelevant and is 
clearly at odds with the infinite sum prediction. 

B. Higgs production 

To evaluate the contribution to the hgg coupling coming from the individual KK fermion 
modes we proceed as in the previous subsection. We now use Eq. (41), and sum up to the 
maximum number of KK modes desired. We perform the sum numerically and show the 
results in Figure 3. Again we are interested in the variation of the couplings with (3 and 
compare them to the result for the Ch gg obtained by calculating the infinite sum, as shown 
in the previous section. 

The two panels of the figure show the contribution to the hgg coupling coming from a 
5D up-like quark (left panel) and the contribution coming from a 5D top-like quark (right 
panel) for (3 values up to 100. We can see how the sums over different maximum number of 
KK modes converge to the infinite sum limit as we vary (3. The approximation obtained by 
considering just a few KK modes is much better for low values of (3. For example, from the 
left panel of Figure 3, for (3 = 2 — > 5, 8 KK modes saturate some 90% of the infinite sum, 
while for (3 = 20, 8 KK modes saturate some 60% of the infinite sum. For (3 = 100 (corre- 
sponding to a Higgs highly localized towards the brane), 8 KK modes represent only some 
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10% of the total KK contribution. This dependence on (3 is in agreement with the results 
found in [16], in which a brane localized Higgs was considered (i.e. (3 = oo) and the first 
few KK fermions considered were found to give a negligible contribution to the hgg coupling. 

We conclude that in all the previous calculations (the top quark Yukawa shift and the 
contributions to the hgg coupling coming from up-type and top-like 5D quarks) we have 
observed the same feature, namely that in the case of a bulk Higgs (small f3), the effect 
of the heavier KK modes decouples (i.e. performing the infinite sum is equivalent to sum 
only over the first few KK modes). On the other hand, when f3 is very large, the heavier 
degrees of freedom do not seem to decouple hinting towards some type of UV sensitivity 
of the brane Higgs case. This is not that surprising since the thickness of a Higgs being 
crushed against the brane is becoming smaller and smaller, and the scale associated with 
the Higgs localization eventually becomes much larger than the cut-off of the scenario. We 
will now see how adding a type of higher derivative operators will be sufficient to make the 
finite sums consistent with the infinite sum results obtained earlier. 

V. THE EFFECT OF HIGHER DERIVATIVE OPERATORS 

We have just seen how the results obtained in the previous section (IV), where we sum 
over a few KK modes agree with the complete KK tower summation of section III only in 
the case of a bulk Higgs boson. When the Higgs is on the brane, or very much pushed 
towards the brane, the results for the two approaches do not seem to agree (see Figure 3 
when (3 — > 100). We will reconcile the two methods by including, in the effective theory 
calculation, the contribution of higher derivative operators. 

In particular we consider the effect of the following operator in the action with a dimen- 
sionful coupling constant Y R (flavor indices are suppressed), 



The operator is of Yukawa-type as it couples two fermions with the Higgs, but it involves 
derivatives of fields. The coupling Y R should be in units of A, the cut-off of the theory, and 
so obviously this operator is cut-off suppressed (we note that the standard 5D Yukawa Y u 
coupling is also dimensionfull and cut-off suppressed, but by two units less than Y R ). Since 




(50) 
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Qr(z) and Ul(z) satisfy Dirichlet boundary conditions on the IR brane, their derivatives 
along the extra dimension can be large after electroweak symmetry breaking and so we focus 
on the operator 

S D J d A xdz (^j [Y R d^ R Hd z u L + h.c] , (51) 

which includes only the wrong chirality fermion components Qr(z) and Ul(z) as it could 
lead to potentially large effects. 

As explained in the previous sections we can proceed in two ways in order to compute 
the effects of this operator. We could study the effect of the operator into the 5D equations 
of motion after electroweak symmetry breaking (ESB) and calculate its effects from these. 
Alternatively, we could solve the equations of motion and perform the dimensional reduction 
before ESB, and then consider the effects produced by the operator by working in this gauge 
eigenbasis. Both methods should be equivalent, but we will follow the second one. In this 
approach, we obtain the effective 4D theory and since it is non-renormalizable, we cut-off 
its spectrum at the cut-off scale thus effectively we only allow a few physical KK modes into 
the calculation. The effects from higher modes are integrated out and encoded in all higher 
order operators of the theory with their effects under control by the cut-off suppression. In 
the case of the Y R operator the potentially large derivatives of Qr(z) and Ul{z) can offset 
the cut-off suppression and so we should keep this operator in the calculations. 

In the approach in which the KK modes are in the gauge basis, the Y R operator will 
affect the fermion mass matrix from eq. (34), and in particular it will contribute to the 
Yjj l q terms. Its effects can therefore be tracked into the effects of these wrong chirality 
terms, as was already noted in the appendix of [20]. We can thus formally treat the situation 
as before, where a truncated version of the infinite mass matrix of eq. (34) is considered (with 
just a few KK levels), but now we redefine the terms Y^ lQr to include the contributions 
from Y R as 

Yv lQr = £ dz (| ) 5 V -M_ (y"U l (z)Q r (z) + Y^ 2 d z U L (z)d z Q R (z^ . (52) 

It is now easy to compute numerically the new effects since from here we just have to 
repeat the previous procedure. The results are shown in Figures 4 and 5. In both figures 
we show the individual contributions coming from the normal Yukawa coupling Y u , from 
the new Y^ coupling, as well as the combined effect. This combined effect is represented 
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FIG. 4. (color online). Contribution to the coupling \cf igg \ (relative to the Standard Model) as a 



function of the Higgs localization parameter f3 when considering only a five-dimensional up-type 
quark, and computed with the higher derivative term discussed in the text in addition to the 
standard 5D Yukawa coupling term. Since both contributions have independent phases we add 
and subtract their generic size to obtain the shaded region of possible values. These results are 
calculated by using only the first 3 KK modes (i.e. considering an effective theory with a cut-off 
of the order the the fourth KK mass). The dimensionless 5D Yukawas are fixed at Y = 2 and the 
KK scale is set at — = 1000 GeV (the overall effect scales as Y 2 v 2 R' 2 ). 

by the shaded region, the reason being that the two types of couplings Y u and Y% have 
independent phases and so can add up constructively or destructively, or in between. In 
Figure 4 we focus on the contribution to hgg due to an up-like 5D quark. In Figure 5 we 
show the predictions for the both shift in the SM top quark Yukawa coupling as well as the 
prediction for the contribution to the hgg coupling coming from a 5D top-like quark. As we 
can see, the shift in the top quark Yukawa coupling can be quite large, and for low values 
of the Higgs localization parameter /3 the shift obtained always results in a suppression in 
the Yukawa coupling. For large values of (3 the shift can be in either direction (suppression 
or enhancement). In the case of the hgg coupling, we see that the contribution represents 
an enhancement with respect to the SM prediction for small values of /3, and again for large 
values of j3 the Ch gg coupling can be either enhanced or suppressed depending on the relative 



18 



5d top-like fermion 



1.03 



5d top-like fermion 



l/R'=1000GeV 
y=2 



1.02 



l/R'=1000GeV 
7=2 



1.00 



0.90 - 



0.95 




0.99 



2 



5 



10 



20 



50 



100 



0.98 ^ 

2 



5 



10 



20 



50 



100 



FIG. 5. (color online). Shift of the top quark Yukawa coupling (left) and contribution to the cou- 
pling Chgg (right), relative to the Standard Model, as a function of the Higgs localization parameter 
/3, when considering only a five-dimensional KK top quark, and including in the computation the 
higher derivative term discussed in the text, in addition to the standard 5D Yukawa coupling term. 
The contributions from each term have independent phases and so we add and subtract their ab- 
solute value to obtain the shaded region of possible values. These results are calculated by using 
only the first 3 KK modes (i.e. considering an effective theory with a cut-off of the order the the 
fourth KK mass). The values of the 5D Yukawa Y u and of Yr are fixed at Y = 2 and the KK scale 
is set at ^7 = 1000 GeV (the overall effect scales as Y 2 v 2 R' 2 ). 

phases between Yr and Y u . For this comparison we have taken the absolute value of both 
couplings and Y u to be the same, i.e. Y = 2, in appropriate units of the cutoff. The main 
feature to remember is that the effects of the higher derivative operator Yr are subdominant 
for small /3 but become dominant for large /3. The large contribution obtained at large f3 
is precisely what makes these new predictions consistent with the results obtained with the 
original infinite sum, and so the higher derivative operators that we have considered here 
somewhat encode the UV sensitivity found in the previous section. 

VI. MISALIGNMENT BETWEEN HIGGS VEV AND HIGGS PROFILE 

In this section we present a discussion on how to treat the case where the Higgs profile 
is different form its VEV profile. This is equivalent to consider the mixing effects between 



19 



the massless zero mode Higgs boson, and the heavy KK Higgs modes and its effects on the 
Higgs observables computed in this paper. 

We follow closely an argument by Azatov [23] and for simplicity we will discuss a simple 
situation in which the 4D effective theory contains only two new heavy vector-like fermions, 
Q and U, doublet and singlet of SU(2)l respectively. This is the situation one would have 
when the KK fermion towers are truncated after the first KK excitation. 

Let's first define our notation for the following quantities 

= Jdz (^J ^ ^ hp(z), (53) 

where v 4 is the SM Higgs VEV and vp(z), hp(z) are the 5D profiles of the Higgs VEV and 
the Higgs physical field, which are generically different. That is, after EWSB, the Higgs 
field is expanded around the nontrivial VEV vp(z) as 

H(x,z) =vp{z) + h(x)hp{z) + .... (54) 

In the case of the bulk Higgs sector considered here, both profiles vp(z),hp{z) are almost 
the same, (see eq. (14)), the order of the misalignment between them being controlled by 
powers of (mhR') 2 (a small quantity). 

We consider all the possible couplings between the Higgs and the fermions of the effective 
theory which after EWSB can be written as the matrix M.{v±, h) as 

[ Y Lu R + o y£ Uh v 4 + Xl v h{x) \(u R ^ 



(ql,QlUl) 



Qr 

\ u »; 



Y Lu R "4 + X* LUH h(x) M Q Yl UR v, + X0 LUH h(x) 

The coupling between the physical Higgs and the two gluons is controlled by the physical 
Yukawa couplings Y? hys and masses Mf hys of the heavier physical fermions running in the 
loop (top quark and KK modes), i.e. 

VPhys 

E jfc* = ^(Y^M-JJ - £ K (56) 

heavy i light 

where Y phys is the physical Yukawa coupling and M phys is the physical fermion mass matrix 
of the setup. Because the trace is invariant under unitary transformations, we can rotate to 
the gauge basis and write 

Tr(Y phys M p i s ) = Tr(Y gauge M-^ ge ), (57) 
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,(55) 



and note that we can now relate this to the matrix M.{v±, h) as 

^(Yg^M-^J = d h log{DetM(v 4 ,h))\ h=0 . (58) 

The procedure is the same as was followed in Section III, i.e. we compute the the determinant 
by expanding in powers of v 2 /M 2 and after combining everything we obtain 

yphys light 

Y,^ = ^ogDetM{v^h)-^ 

heavy » 

I ^QlUrUlQr _ ^UlQrQlUr \ /rg\ 

4 \ m q Mu m q Mu y K ] 

This result is the equivalent to eq. (40) with the effect of the misalignment between vp(z) 
and hp(z). One sees that the difference lies in the substitution of one of the Y terms by an 
X term, and so the correction to the result of eq. (40) is 

/ \ X QlUr ~~ Y Q l Ur ) Y U l Qr \ X U l Qr ~ Y UlQr) Y QlUr\ 

Sc "» = v ' ( l^h, JVA, J ' (60) 

which is controlled by 

{ X Lqr- Y $ lQ r) =j dz (f) U L {z)Q R {z) (hp{z) - ^) , (61) 

and 

{ X Lur - Y Lu R ) =jdz (I)' 'Ql(z)U r (z) (hp{z) - ^) , (62) 
and since the misalignment between vp(z) and hp(z) can be computed perturbatively [20] 



as 



h " {z) -^- — \WW) + WW) + { h )l (63) 



vp(z) vp(z) ( m 2 h R- mi : 

npyt) - 

we obtain 



/ Y p Y p+2 \ 

(x uo -Y uo )= m 2 h R' 2 4t^- + - UlQr + 0(miR' 4 ). (64) 

V u lQr u lQr J h ^ 2 (4 + /3) 4v/(l + /3)(3 + /3) J y h > y > 

In other words, the effect of considering the misalignment between vp(z) and hp(z) is to 
add a correction with the same structure as the result of eq. (40), but with a suppression 
of (m h R) 2 , i.e. the correction is at most (9(1%), and becomes much smaller for increasing 
values of 8. 5 



5 The dependence on f3 of the integrals ^ij l q r and Yq lUr is quite mild and so, in terms of order of magnitude, 
wehave< Qfl ~<+ g 2 fi . 
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VII. DISCUSSION AND OUTLOOK 



In this work we have presented the results for the predictions of Higgs phenomenology in 
a toy-model RS setup in which the Higgs field is allowed to propagate in the bulk and with 
a single 5D fermion field. Our results can be extended to three families to include full flavor 
effects, but the generic predictions that we would obtain are expected to be basically the 
same as the ones presented in [18, 20]. That is, that in the context of flavor anarchy, where 
the action parameters are all of the same order but with more or less random values and 
phases (with the constraint of obtaining correct SM predictions) the couplings of the Higgs 
with fermions and gluons and photons can receive important corrections, either enhancing or 
suppressing the SM predictions. However, the two references mentioned present calculations 
performed by including the effect of all the KK fermions, technically assuming an infinite 
cut-off for the model (where a brane Higgs is considered). In general, all these scenarios 
break down at a low cut-off, becoming strongly coupled for both gauge and Yukawa inter- 
actions. The implicit assumption made in [18, 20] was that the effects of the heavier modes 
should decouple quickly, at least for the case of a bulk Higgs field. The main motivation 
to perform the calculations by considering the full infinite fermion KK tower, as well as 
pushing the Higgs into the brane was mainly of technical nature. Indeed both the flavor 
structure of the Higgs Yukawa couplings as well as the coupling to gluons and photons can 
be computed analytically with those ingredients. In [20], the authors checked analytically 
that the corrections to the Higgs Yukawa couplings were actually of the same order for a 
bulk Higgs and a brane Higgs. 

However it was pointed out in [17] that in the brane Higgs case, the effects of the heavier 
KK fermion modes do not decouple and that they all contribute evenly in the computation 
of the Higgs couplings in the model. On the other hand, we showed in sections III and 
IV of this paper that the heavier KK modes in the case of a bulk Higgs do decouple very 
quickly, so that the analytical result obtained by using the infinite KK tower approaches 
with great precision the numerical result obtained by considering an effective theory with 
only a few KK fermion modes. Moreover, when considering the effective theory with only 
a few KK modes, one should include in the action all possible operators and in particular 
the higher derivative ones introduced in Section V. These effects were omitted in [17], and 
as we showed in this work, the importance of these operators increases as the Higgs is more 
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and more localized towards the brane. In [16], the authors considered an RS setup with a 
highly localized Higgs and the presence of only a few KK fermions and studied the effects 
on the Higgs couplings to gluons and photons, among other observables. In the limit of 
the SM gauge group (they did consider an extended gauge group) they found no significant 
deviations from the SM predictions. Indeed this result is consistent with our findings of 
Section IV (no higher derivative operators invoked yet), since as it can be seen on Figures 5 
and 4, the shift in Higgs Yukawa couplings and the new effects to Higgs-gluon-gluon coupling 
vanish in the limit of highly localized Higgs (large (5 parameter). On the other hand, in [17] 
it is claimed that large effects should be present in the case of a brane Higgs and with only a 
few KK modes present in the effective theory (and no higher derivative operators), a result 
inconsistent with both our findings and those found in [16]. We can trace the origin of 
the disagreement in their calculation of the Higgs Yukawa couplings. Those are computed 
by using the full 5D equations of motion, which as we have said earlier is equivalent to 
considering the complete tower of KK modes. Then, using these couplings, they calculate 
the hgg radiative coupling but now including only a finite amount of KK fermions. This 
treatment leads to a highly suppressed top quark Yukawa coupling (due to effects from 
the infinite KK tower) and a vanishing contribution to hgg from the loops of KK fermions 
considered (one would need the whole tower to obtain a finite effect). Their end result is a 
suppressed top quark Yukawa coupling and a suppressed hgg coupling (due to the smaller 
top quark Yukawa), predictions which are at odds with the findings of [16, 18] and of this 
paper. 

The procedure of [17] seems inconsistent because essentially the authors use infinite KK 
degrees of freedom in one part of their calculation (the SM quark Yukawa couplings compu- 
tation via equations of motion) but then they truncate the KK degrees of freedom in order 
to compute the hgg coupling. In any case, had they included the higher derivative operators 
introduced in this paper, their results would have changed dramatically since then, the effect 
to hgg coming from the top quark Yukawa loop would remain basically the same, but the 
effects due to loops of a few KK fermions would dominate the overall effect (and thus the 
result would start to become consistent with the findings of [18]). 

Also, the predictions of [16] should change if one considers the effects of the higher 
derivative operators introduced in Section V. In that situation, the Higgs couplings can 
receive large corrections, and can be of any sign (suppression or enhancement) due to the 
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different phases present in the couplings Y and Yr. In fact we have found here that for a 
Higgs field in the bulk, our results are more predictive than for a brane Higgs field, because 
the effect of the higher derivativer operators is subdominant for a bulk Higgs field. 6 The 
effects from only the 5D Yukawa operators are aligned [20], and thus all the KK quarks 
add up in phase. In that situation we can have definite predictions for the effects caused 
by a single family of fermions, i.e. it will produce a suppression in the light quark Yukawa 
coupling and an enhancement in Higgs boson production (as well as suppression in the 
Higgs to photons coupling) [18], with the caveat of taking the dimensionless couplings of 
both Yukawa terms and higher derivative operators to be the same (consistent with the usual 
assumption that all 5D coefficients have to be of the same order). Taking into account the 
three fermion families in conjunction with a bulk Higgs field might weaken this prediction 
due to complicated flavor mixings and structure, but still one should be able to draw a 
correlation between Yukawa couplings and Higgs production (and h — > 77) for the case of a 
bulk Higgs field. The parameter space of the bulk Higgs scenario can therefore be under a 
tighter pressure as more and more precise experimental measurements in Higgs observables 
at the LHC become available. In particular if the predicted and correlated deviations of 
Higgs couplings is not clearly observed this should put bounds on the KK scale of the bulk 
Higgs scenario. 

The situation for a Higgs on the brane is different. The higher order derivative terms 
are now important. Each KK tower of light quarks and the top will contribute to the 
hgg coupling, but their effect depends on arbitrary relative phases (between Y R and Y 5d ), 
and so one cannot make a firm statement about the magnitude and phase of the overall 
contribution: it can be a suppression or an enhancement, or in between. 

Finally we comment again on the apparent problem of a highly localized Higgs scenario 
(brane Higgs) in which predictions made from a truncated fermion KK tower are very differ- 
ent from predictions made from an infinite fermion KK tower. This apparent UV-sensitivity 
can actually be lifted by considering the higher derivative operators described here (first in- 
troduced in [20]). When these are included, the predictions made with a finite KK fermion 
tower become consistent with the original predictions obtained with an infinite fermion 
gain, the reason for this is that the value of the derivatives of the bulk fermions is suppressed by the higher 
value of the 5D cutoff. When the Higgs boson is pushed towards the brane, the derivatives of these fermions 
fields (with the "wrong" chirality) becomes larger and larger, and the 5D cutoff does not suppress anymore 
the effect of these operators. 
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tower. A more esthetic problem with the brane Higgs scenario remains, since the definition 
of the Higgs operators seems highly unnatural, if one understands a brane Higgs field as a 
limit of a bulk Higgs field. All operators involving Higgs fields will have to have a precise 
and definite dependence on (3 (a large number), which seems quite contrived, specially in a 
framewrok in which no big numerical hierarchies should arise from fundamental 5D coeffi- 
cients. In any case, with the ansatz outlined in the text and reviewed in Appendix B, one 
can still work consistently with a brane Higgs field as a limit case of a bulk Higgs field. 
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VIII. APPENDIX 



Appendix A: Some Explicit Analytic Results 



A 1 Y 

From equation (39) the shift defined as = , can be also derived from 

mv 4 t> 4 m 



yti yu* yam 

^4 \^ I qU j I U j Q i I Q i u / a 1 \ 



™± Y- *rf Mo Mr 



1,3 



Therefore simply replacing the second Yukawa coupling with the Yukawa coupling of the 
Higher derivative operator, Y R will give 

Here, we present explicit analytic expressions for the hgg production and also the Yukawa 
coupling-mass shifts by performing the infinite sums over the KK modes. We also include 
the result given in reference [20] for the shift due to the usual Yukawa term, Y d QHU, for 
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completeness. 7 To summarize, we have 

■(l- e l-2c,)(l- e l+2*.) 



c hgg 



+ 



2m d Rl 2 2 + c u -c q + /3 
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for the hgg production and 
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for the shifted Yukawa coupling. Also, there is a misalignment due to the kinetic term [20], 
which as discussed in the text, is only important for the case of the third generation quarks. 
We do not repeat that result here. For the higher derivative term the shift is: 



(A 
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Appendix B: Prom Bulk to Brane 



We summarize the matching prescription for operators containing Higgs field for the case 
where the Higgs boson is localized on the brane. As explained in Section II, these pre- 
scriptions insures that the 5D bulk Higgs scenario transitions smoothly to a brane-localized 
Higgs case. The brane prescription for the Higgs associates a delta function to the Higgs 
normalization integral 



7 We have reproduced this result using the equation (A2), and our results match the one given in the text of 
[20]. Note however that there a few typos in the eq. Al of their appendix. 
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f(%[M^ = i- 

Jr z 

As the HH, rather than H field, is associated with a 5 function, one must include a (3 
dependence to the bulk Higgs fields to be able to match operators, in the limit /3 — > oo to 
the brane ones. The conversion is: 

HH — >■ HH , 
#M ->■ -^=HHH , 

for matching brane to bulk in the appropriate limit. 

So for the shift, we have contributions from Y 2 and Yr. As we are dealing with an effective 
theory, we look at the effect of summing over a finite number of modes, let's say 3 to 5. 

For the case of brane Higgs, the contributions for a finite number of modes for Y 2 give 
exactly (because of boundary values on the brane). This confirms the work of [16]. How- 
ever, we must add higher order operators Yr, which give a significant result (converging to a 
constant for f3 — > 1000 and anything beyond). The result obtained by summing over a finite 
number of modes in the brane on the Yr contribution must be compared with the result in 
the paper by [20] for the infinite sum of Y 2 on the brane. 

For bulk Higgs, the shift contribution from a finite number of modes on the Y 2 contri- 
bution is no longer 0. However, adding to this the Yr contribution, we notice that the Yr 
contribution for bulk Higgs is much smaller (two orders of magnitude) than the correspond- 
ing one in the brane. This is a clear indication that higher order corrections are much more 
important for the brane Higgs case than for the bulk. 
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